g-Breathers and thermalization in acoustic chains with arbitrary nonlinearity index 
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Nonlinearity shapes lattice dynamics affecting vibrational spectrum, transport and thermalization 
phenomena. Beside breathers and solitons one finds the third fundamental class of nonlinear modes 
- g-breathers - periodic orbits in nonlinear lattices, exponentially localized in the reciprocal mode 
space. To date, the studies of g-breathers have been confined to the cubic and quartic nonlinearity 
in the interaction potential. In this paper we study the case of arbitrary nonlinearity index 7 in 
an acoustic chain. We uncover qualitative difference in the scaling of derealization and stability 
thresholds of g-breathers with the system size: there exists a critical index 7* = 6, below which both 
thresholds (in nonlinearity strength) tend to zero, and diverge when above. We also demonstrate 
that this critical index value is decisive for the presence or absense of thermalization. For a generic 
interaction potential the mode space localized dynamics is determined only by the three lowest order 
nonlinear terms in the power series expansion. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 
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Nonlinearity is a key player in a number of fundamen- 
tal dynamical and statistical physical phenomena includ- 
ing thermal conductivity, wave excitation and propaga- 
tion, electron and phonon scattering. Its impact is of- 
ten counterintuitive. E. Fermi, J. Pasta and S. Ulam 
(FPU) hypothesized that it underlies thermalization in 
crystals and ran the celebrated computational experi- 
ment with oscillatory chains with nonlinearity in coupling 
[1]. Paradoxically, the low wave number q excitations did 
not spread over the spectrum, staying localized in several 
neighbor modes. More than that, it showed almost ex- 
act recurrencies of the energy to the initial mode at large 
times. 

At present we know that this scenario is generic and 
holds in a certain parameter domain. There exist weak 
and strong stochasticity (chaos) thresholds in the en- 
ergy/nonlinearity strength: above the former oscillations 
become chaotic, but stay localized in the mode space; 
only above the latter the oscillations delocalize rapidly 
and thermalization occurs [2-4]. 

Despite certain achievements in the continuum-limit 
approximation with integrable soliton dynamics Q and 
analytic and numeric estimates of the strong chaos 
threshold 

d-El, only the recently developed theory of 
q-breathers (QB) has explained all main ingredients of 
the FPU problem (l2| . 

QBs are exact periodic solutions, continuations of the 
linear modes in the nonlinear regime, exponentially local- 
ized in the (/-space and stable in the suitable parameter 
range. FPU trajectories were shown to be small per- 
turbations near stable QBs, hence the absence of ther- 
malization and recurrencies. The instability of QBs was 
found to coincide with the weak stochasticity threshold, 
and derealization - with the strong stochasticity one and 
thermalization. 

However, there is a principal question that has never 
been addressed: what is the impact of arbitrary nonlin- 
earity in the interaction potential? Indeed, all studies 



were done with the cubic and quartic ones. A generic 
nonlinear interaction potential yields higher order non- 
linear terms in the Taylor series expansion about equi- 
librium. As the magnitude of these terms increases with 
the energy it could be expected that these terms, and not 
the qubic and quartic, determine the derealization and 
stability of QBs, and, broadly, thermalization properties 
of the systems. 

In this Letter we present the theory of QBs in acous- 
tic chains with arbitrary nonlinearity index in the in- 
teraction potential, exemplifying in the FPU chain. We 
demonstrate that there exists a critical nonlinearity in- 
dex 7* = 6, such that for lower order nonlinearities the 
derealization and stability thresholds of QBs approach 
zero as the chain length scales to infinity, while for higher 
orders they diverge. We find that the critical index value 
plays similar role in thermalization of the chain, which 
is rapid below it and non-observable above. It follows 
that only three lowest order nonlinear interaction terms 
are decisive for localization of energy or thermalization 
of generic acoustic chains. 

The FPU system models an atomic chain by N equal 
mass particles coupled by springs with linear and nonlin- 
ear interaction terms: 
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where x n is the deviation of the n-th particle from the 
equilibrium, and the boundary conditions are fixed xq = 
xn + i = 0. The classic a— and ft— FPU models arise for 
7 = 3, 7 = 4 respectively. 

Canonical transformation x n (t) = 

defines 



N 



the reciprocal 



space of N normal modes with the amplitudes standing 
for coordinates Q q (t). Dynamical equations in this case 
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read 
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where w g = 2 sin (■Kq/2(N + 1)) are the normal mode fre- 
quencies, and 
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are the intermode coupling coefficients, defining the long- 
range coupling between g-oscillators. 

In order to establish the continuation of QB in the non- 
linear regime we set x = 0, the energy E qa — Em. the 
normal mode q = qo, while the other q-oscillators are at 
rest. These initial conditions correspond to the periodic 
trajectory in the phase space of the linear system. Fol- 
lowing [13j|, the absence of resonances noj qo ^ ui q ^ qa is 
a sufficient condition for the continuation and the latter 
holds for generic finite N fl4| . 

Following this argument we develop the perturba- 
tion theory for (J^J) in powers of the small parameter 
a = X /[2(N + I)]?/*-*: Q q (t) = ££oCr*Q«(t). The 
frequency of the nonlinear mode is sought in the form 



(t) by means of the Lindstedt-Poincare 



method. Take the low-frequency qo N linear mode 
Qq (t) = Ao cosui qo t, Q q (t) — 0, q ^ qo and its linear 
frequency LJqg = oj qo as a zero-order approximation. The 
structure of the coupling coefficients determines that the 
only non-zero first order terms are Q/^_i)go^ 5 wnere 
k = 4, 6, . . . , 7 for even and n = 3, 5, . . . , 7 for odd 7. 
The coupling coefficients wrap up as -B( re _i) go , go ,... >go ~ 

a 

form of linear forced oscillators: 
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and, finally, the equations of motion get the 
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The amplitude of the response on the resonant frequency 
(k — l)u! qo reads 
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FIG. 1: Localization of QBs in dependence on the chain size 
for the nonlinearity index below and above the critical 7 = 6. 
Here q = 1, X = 0.001, 7 = 4 and 7 = 10 



and the frequency of the periodic orbit is given by 
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Rewriting ([S]) in terms of linear mode energies E q 
w q A q /2, we obtain localization factors A 7iK : 
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One can also define the localization length in the mode 
space £ 7 , K such that £( K _i) 9o = e ( - K ~^ qo ^^- K E qo . Then 
it follows £ 7 , K = — (k — 2)g /(21nA 7jK ). 

Delocalization criterion is, thus, max{A 7iK } ~ 1. It 

corresponds to strong interaction between the linear 
modes and the breakup of perturbation theory. 

The crucial property of (JT]) is the scaling of localiza- 
tion factor A with the system size. For 7 < 6 there is 
a delocalization threshold in N above which the corre- 
sponding QB delocalizes. For 7 = 6 the localization is 
independent of the system size in the leading order ap- 
proximation, while for 7 > 6 the localization strengths in 
longer chains. 

Existence of the critical 7* = 6 is of principal 
importance for the general problem of thermalization 
and equipartition. As indicated above, the thermaliza- 
tion/strong chaos threshold in a- and /3-FPU is directly 
related to delocalization of QBs. It follows, that above 
the critical nonlinearity index the onset of strong chaos is 
retarded to much higher energies; moreover, the thresh- 
old in energy will grow with the system size. The second 
consequence is that for a generic nonlinear interaction 
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FIG. 2: Evolution of participation number in dependence on 
the nonlinearity index. Here the mode go = 1 is excited to 
E = 1 at t = 0, x = 1) N = !000 



potential the thermalization threshold will be fully deter- 
mined by the three lowest-order nonlinear terms (7 < 6) 
in the Taylor expansion about an equilibrium. 

Let us illustrate theoretical predictions with numerical 
results. First, we compute the QBs in ([T]) for different 
nonlinearity index, following the numerical scheme devel- 
oped in [12j. It is clearly seen that localization weakens 
for 7 < 6 and improves for 7 > 6 with the increase of the 
system size N (FigQ]). 

To probe the dependence of the thermalization pro- 
cess on the nonlinearity index directly we perform the 
long-time integration of |T]) with the initial conditions 
corresponding to the lowest mode excited to the energy 
E = 1. The other parameters are N = 1000, x — 1 
and the integration time t — 10 7 is of the order of 5000 
oscillation periods of the seed mode. To characterize 
the number of effectively excited modes we calculate the 
time dependence of the so-called participation number: 
P = l/(X^Li e q)j where e q = E q /E are the normalized 
modal energies. (It is easy to see, that if, say, to modes 
are excited to the same energy and the energy of the oth- 
ers is then P = to) . This measure allows us to visualize 
the mode excitation process and the results clearly show 
thermalization for 7 < 6 and its absence (on the observed 
timescale) for the higher order nonlinearities (FigJ5J). 

Let us study the linear stability of QBs for even non- 
linearity indexes (note that instabilities are suppressed 
for odd 7 = 3 in the whole region of existence of the 
localized QB solution We keep in mind that it is 

tantamount to the onset of weak chaos, when chaotic tra- 
jectories stay localized in the mode space on the compu- 
tationally accessible times, but are believed to thermalize 
through Arnold diffusion, eventually 0, [Hj]. For that we 
linearize equations of motion ^ around the QB solu- 
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For 7 = 4 the strongest instability comes from primary 
parametric resonance in ((5J and involves pairs of resonant 
modes q,p — go ± to, m = 1, 2, 3, . . .. For higher indexes 
the number of potential instabilities increases and one 
needs to take into account all pairs of the type q,p = 
(k — 2)qo ± m, m = 1, 2, 3, ... , K = 4, 6, . . . ,7 — 2. For 
each pair one excludes non-resonant terms and obtains 
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Standard stability analysis [12| reveals that the first in- 
stability to happen corresponds to n = 4 and to = 1 so 
that the instability condition reads: 
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Remarkably, the same critical nonlinearity index 7* = 

6 characterizes the scaling of the instability threshold x* 
with the system size N, as it does for localization. For 

7 < 6 the instability threshold gradually diminishes as 
TV grows; for any fixed set of the energy and nonlinear- 
ity parameters there exist such N above which all QBs 
will become unstable. On the contrary, for 7 > 6 the 
instability threshold grows with the system size; even if 
QBs were unstable in small systems they will attain sta- 
bility if the chain length exceeds a threshold determined 
by (fTTj) . Recalling the relation between instability and 
the onset of weak chaos, we come to similar conclusions 
as to the strong chaos threshold. Namely, for the FPU 
systems with the nonlinearity index higher than the crit- 
ical one, weak chaos gets suppressed in large chains. In 
case of generic nonlinear coupling potential only three 
lowest order nonlinear terms in the expansion about the 
equilibrium influence and cause the development of weak 
chaos regime. 

Finally, we report numerical results on the QB stabil- 
ity. We compute Floquet multipliers for QBs centered at 
go = 3 in dependence on the nonlinearity coefficient tak- 
ing E = 1, N = 16, 32, 64, and 7 = 4, 6, 8. In accord with 
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FIG. 3: Dependence of the absolute values of the Floquet 
multipliers fj, of QBs on the nonlinearity strength for (a) 7 = 4 
and (b) 7 = 8. Circles mark N = 16, squares - N — 32, 
triangles - TV = 64 



tween delocalization/instability and strong/weak chaos 
thresholds, the latter must demonstrate the same regu- 
larities. Indeed, the long-time evolution from the low- 
est mode excitation show fast thermalization below the 
critical nonlinearity index and its absence above. For a 
generic interaction potential, localization of energy and 
thermalization processes are determined only by three 
lowest order nonlinear terms in its power series expan- 
sion; if they are absent, then thermalization is strongly 
suppressed in large systems and extremely high energies 
are required for it to be observed. 

Beside the fundamental physical problems of thermal- 
ization and ergodicity, the results build on the theory of 
nonlinear oscillations and energy localization in spatially 
discrete systems with the high demand coming from ex- 
periments and applications with micro- and nano-scale 
electro-mechanical oscillator arrays lil 17|. 

I thank S. Flach, O. I. Kanakov and K. G. Mishagin 
for useful discussions. 
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FIG. 4: Theoretical solid lines, and numerical, markers, 
instability thresholds in dependence on the nonlinearity index 
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pip, the instability threshold max{/i} > 1 increases with 
N for 7 = 4 and decreases for 7 = 8 (FigOJ) . The quanti- 
tative agreement between the theoretical predictions for 
X* and numerics is very good (Fig(4|). 

In conclusion, we develop the theory of q-breathers 
in acoustic chains with arbitrary nonlinearity index in 
the coupling potential. Theoretical and numerical stud- 
ies of the mode space localization and stability reveal 
that there exists the critical nonlinearity index 7* = 6, 
below which q-breathers lose localization and stability 
as the system size is increased. In a marked contrast, 
these properties strengthen in larger chains for higher 
order nonlinearities. According to the direct relation be- 
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